Nodal quasiparticles and the onset of spin density wave order in the cuprates 
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We present a theory for the onset of spin density wave order in the superconducting ground state 
of the cuprates. We compute the scaling dimensions of allowed perturbations of a 'relativistic' fixed 
point with 0(4)® 0(3) symmetry, including those associated with the fermionic nodal Bogoliubov 
quasiparticles. Analyses of up to six loops show that all perturbations with square lattice symmetry 
are likely irrelevant. We demonstrate that the fermion spectral functions are primarily damped 
by the coupling to fluctuations of a composite field with Ising nematic order. We also discuss the 
influence of quenched disorder. 



A large number of experiments have investigated the 
magnetic correlations in the cuprate superconductors [l[ . 
There is good evidence for collinear spin density wave 
(SDW) correlations near the wavevectors 
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where a is the square lattice spacing and the incommen- 
surate "d varies with the hole doping concentration, 8. 
Especially in the La based compounds, this order is best 
formed near 8 = 1/8, where d ss 1/8. This SDW order 
coexists with superconductivity (SC), and a theory 0] for 
the tuning of the transition between the SDW+SC and 
SC states by an applied magnetic field has been tested 
in a number of experiments [1, 0> S, 0| ■ 

In this paper, we will examine the interplay between 
the SDW order and the fermionic Bogoliubov quasipar- 
ticle excitations of the superconductor. Because of the 
d-wave nature of the Cooper pairs, these quasiparticles 
have the spectrum of massless Dirac particles whose en- 
ergy vanishes at four wavevectors (the nodal points) in 
the Brillouin zone (±Q,±Q); these nodal wavevectors 
bear no relation to the SDW ordering wavevectors K1.2. 
Using a theory for the critical fluctuations of the SDW 
order developed in Refs. @, H, 0], we will show that the 
dominant coupling between nodal quasiparticles and the 
SDW order is to a composite Ising nematic order [loj ]. 
associated with the choice of the SDW ordering between 
the Ki and K2 wavevectors. We emphasize that the 
unique selection of Ising nematic order here is not a phe- 
nomenological assumption, but a non-trivial consequence 
of the internal structure of the SDW fixed point of De 
Prato et al. @; indeed, this is the main point of this 
paper. Our theory predicts nodal quasiparticle spectral 
functions which have been measured in numerous pho- 
toemission and scanning tunnelling microscopy experi- 
ments. The importance of the Ising nematic order in 
damping nodal quasiparticles was noted recently by Kim 
et al. [11] in a different model [T3 |. and we will compare 
our results to theirs. 

Our theory is expressed in terms of the complex- valued 
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FIG. 1: Transformations of the fields under operations which 
generate the symmetry group: T XtV = translation by a lattice 
spacing in the x, y directions, R = rotation about a lattice site 
by 90°, I — reflection about the y axis on a lattice site, and 
T = time reversal. The theory is also invariant under spin 
rotations, with i a vector index and a, f3 spinor indices. We 
define T as an invariance of the imaginary time path integral, 
in which transform as the complex conjugates of $i,2i, 

while ^1 2a are viewed as independent complex Grassman 
fields which transform as 2a — > &\ 2c tV ■ 



order parameters $u and $2i, where i — x,y,z denotes 
spin components, which are related to the spin operator, 
Si, at position r and imaginary time t by 

Si(r, t) = Re [e* K - r $ H (r, r) + e* K - r $ 2l (r, r)] . (1) 

This parameterization implies the transformation prop- 
erties of $1,21 under the symmetry operations on the un- 
derlying square lattice Hamiltonian, which are summa- 
rized in Fig[TJ By writing down all terms invariant under 
these operations, and expanding in powers and gradients 
of $1,21, we obtain the following quantum field theory 
for the fluctuations of the SDW order 0, @] with action 
J cPrdTC® and 
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where are spin-wave velocities, r is the coupling 
which tunes the system across the transition to a state 
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with SDW order, and 1*1,2 and 101,2,3 are the crucial quar- 
tic couplings which stabilize and select the SDW order. 
Linear spatial derivative terms such as Q\d x Q\ are also 
permitted, but can be absorbed by a redefinition of the 
incommensuration d. For some of the cuprates with 
5 ~ 1/8, the ordering is commensurate with ■& = 1/8; 
in this case a lock-in term ~ (^f) 4 + (^l) 4 is allowed, 
but this eighth order term is clearly irrelevant near the 
critical point, and can be safely neglected in our consid- 
erations here. 

The weight of the experimental evidence, reviewed e.g. 
in Ref. [l], is that the SDW order in the cuprates is 
collinear in spin space, and picks a definite spatial di- 
rection by a choice of condensing either <!>i or $2. So a 
particular state has 
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or an equivalent state with 1 <-> 2, with m an arbitrary 
real vector and e l8 a common complex phase for the com- 
ponents of <&. In mean field theory, such a state is the 
ground state of the theory <S$ in a subset of the region of 
parameters satisfying r < 0, w\ > 0, w± + W2 + W3 > 0, 
U2 < Min[0, to 1 — ui,wi + u>2 + W3 — ui,(wi + w 2 — 
ui)/2, (wi + W3 — ui)/2]. Going beyond mean field the- 
ory, a detailed analysis of the transition into this SDW 
state under the Lagrangian was provided in Ref. 0. 
It was found that such a transition was in the domain 
of attraction of a fixed point with 0(4)<g>0(3) symmetry 
within two three-dimensional perturbative schemes: (?) 
the massive zero-momentum scheme (MZM) defined in 
the unbroken phase (to six loops), and (ii) the minimal 
subtraction scheme (MS) without expansion in 3 — d (d 
is the spatial dimension, and d = 2 was set after renor- 
malization) defined in the massless critical theory (to five 
loops). High-order perturbative calculations indicate the 
stability of this point also in the SDW model i.e. 
the quartic perturbations present in £$ turn out to be 
irrelevant at 0(4) ® 0(3) fixed point. The values of the 
quartic couplings in the fixed-point Lagrangian L*^ obey 
w* = u* — 1*2 and w 2 = =11%, and the spin- wave 
velocities are equal v\ = The symmetry of the fixed 
point becomes explicit by introducing a 12-component 
real field (p a i, with a = 1 ... 4, related to the SDW order 
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now the 0(4) and 0(3) rotations act on the a and i 
indices respectively. The critical exponents associated 
with this fixed point are the correlation length exponent, 
v = 0.9(2), and the anomalous dimension of the SDW 
order parameter, 77 = 0.15(10). 

We are now ready to turn to our new results on the 
coupling between the above SDW fluctuations and the 
excitations of the superconductor. 

Density fluctuations of the superconducting conden- 
sate will couple to l^il 2 + I <E>2 1 2 - For short-range repul- 
sive interactions between the Cooper pairs, this coupling 



has scaling dimension [13[ (2 — 3^)/2 < and so is ir- 
relevant. The long-range Coulomb interactions further 
suppress density fluctuations, and so we need not con- 
sider this coupling further. 

Far more interesting, and subtle, are the couplings be- 
tween $12 and the fermionic nodal quasiparticles. Let 
us denote the electron annihilation operator with mo- 
menta in the vicinity of the nodes as (Q,Q), (—Q,Q), 
(-Q,-Q), and (Q,-Q) by f la , f 2a , f 3a , and f 4a re- 
spectively; here a =1,1 is an electron spin index. Next, 
we introduce the 4-component Nambu spinors 
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where e a p = —ep a and £fj_ = 1. We will use Pauli ma- 
trices T l which act on the Nambu particle-hole space, 
while <j l a p will act on spin space. In the vicinity of the 
nodal points, we can expand the standard Bogoliubov 
Hamiltonian of a d-wave BCS superconductor in term 
of gradients of ^1.2, and obtain the low energy theory 
Sqi = J d 2 rdr£q, with 

= *I (d T - i^={d x + d y )r z - i^(~d x + d y )rA *! 



* 2 ( d T - i^(-0 X + dy)T Z - i^(d X + 8y)T X ) * 2 . (6) 
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Here vf and va are the fermionic velocities normal and 
parallel to the underlying Fermi surface. As was the case 
with 5$, the transformations of the fermionic fields ^1.2 
are also crucial to our analysis, and these are summarized 
in Fig. [IJ 

We now write down the most general couplings be- 
tween the SDW fields $1,2 and the Dirac fermion excita- 
tions $1,2 which are allowed by the symmetries in Fig. [1] 
The translational symmetries T XtV immediately rule out 
any terms linear in i>: the SDW fluctuations scatter the 
electrons by wavevectors K 1,2 and these do not, in gen- 
eral, connect the low energy excitations at any pair of 
nodal points Moving on to terms quadratic in $, 

which are generated by a virtual scattering process to 
a high energy intermediate fermion excitation (see also 
Ref. [l4j), all the allowed terms are 



£i = Ai(|$i| 2 
£ 2 = A 2 (|$i| 2 

£3 = tijk 
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Key to our remaining analysis are the scaling dimen- 
sions of the A couplings at the 0(4)<g)0(3) fixed point C% 
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discussed earlier, at which the fermions are decoupled 
from the SDW order. The free Dirac fermion theory 
is also invariant under the same scaling transformation 
as with dim^] = 1, and this allows us to obtain the 
dim [A] from a knowledge of the dimensions of the associ- 
ated composite operators of 4> under In particular, 
we perturb by the same quadratic $ operators above 
but without the fermion terms, and subtract 2 from the 
scaling dimension of such a. perturbation. So e.g. we 
consider the theory £J + Ai (|$i| 2 + |$2| 2 ) and obtain 
dim[Ai] = dim[Ai] — 2. We classify such perturbations 
under representations of 0(4)<g>0(3), with Dr- indicating 
a dimension D representation of O(K). Then Ai corre- 
sponds to a (14,13) operator, A2 to (94,13), and A3 to 
(64,33). 

The (I4, 13) operator tunes away from the transition 
by changing r, and so 



dim[Ai] = l/v -2. 



(8) 



From the value of v quoted earlier, we conclude that Ai 
is safely irrelevant. 

The second coupling between $ and the fermions in- 



volves the operator 



$,12 



$2! i which is the Ising 



nematic order parameter 1(J> 15 1 ■ This measures fluctu- 
ations which break the orthorhombic symmetry, C^ Vl of 
the square lattice to tetragonal symmetry, C^, by mak- 
ing a choice of the SDW order between the Ki and K 2 
wavevectors. This operator is a component of the (94, 13) 
operators, and to obtain its scaling dimension we com- 
puted and analyzed perturbative series within the MZM 
(to six loops) and MS (to five loops) expansion schemes. 
Details on the perturbative calculations and the series 
will be reported elsewhere. The comparison of the anal- 
yses within the MZM and d = 3 MS schemes represents a 
non trivial check of the results. In order to estimate the 
scaling dimension, we need to resum the series and then 
to evaluate them at the fixed point values of the quar- 
tic couplings. The resummation of the series is performed 
by using the so-called conformal mapping method, which 
exploits the knowledege of the large-order behavior of the 
expansions 1^, 17 1. In this manner we obtained 



dim[A2] = 



1.95(18) - 2 
1.90(27) - 2 



MZM 
rf = 3MS 



(9) 



These results favor the irrelevance of A 2 at the decoupled 
fixed point, although the precision of our results does not 
allow us to state this conclusively. In any case, it is clear 
that |dim[A 2 ]| is nearly zero, indicating that A 2 hardly 
flows under the renormalization group (RG). 

Finally, the A3 , A 3 couplings in Eq. ([7|) involve compos- 
ite $ operators (which are components of (64, 33)), which 
measure spiral spin correlations, as is easily deduced from 



Eq. |T]). By an analysis as above, we found 



dim^A^/ 1 - 16 ; 8 ;- 2 

1 31 \ 1.24(8) - 2 



MZM 
d = 3MS 
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indicating that such perturbations are clearly irrelevant. 

We have now established one of the main results of our 
paper: the only important coupling of the nodal fermions 
to the SDW order is A 2 , which couples the fermionic and 
SDW contributions to the Ising nematic order <fi. 

Let us now compute of the influence of the SDW fluctu- 
ations on the fermionic spectral functions perturbatively 
in A2. At leading order, we use the propagator for 0, 
implied by its scaling dimension in Eq. ^ 



G<£(p,w) 



1 -l/2-dim[A 2 ] 
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where p and ui are the momentum and imaginary fre- 
quency carried by 0. To second order in A 2 , for the 
fermion Green's function G$ we obtain the self energy 
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From the finite temperature (T) generalization of 
Eq. (fl"2"|) , we obtain at the nodal point, the retarded, 
real frequency self energy lm£*(0,0) ~ yi-2dim[A 2 ] in 
the quantum-critical region fl2| . This is the leading T 
dependence for dim[A2] < 0. For dim[A 2 ] > 0, the RG 
flow must be computed to higher order in A 2 , and if there 
is a fixed point with a non-zero A 2 , then T is character- 
istic energy scale for quasiparticle damping. Given the 
near-marginality of A2 (Eq. ©), it is valid to ignore its 
RG flow in experimental applications, and so we obtain 
a quasiparticle scattering rate which is practically linear 
in T . Similarly, at the nodal wavevectors, a nearly linear 
frequency dependence is found in the fermion self energy 
at T — 0. These are phenomenologically attractive fea- 
tures [m. 

In describing the q dependence of the fermion spec- 
trum, we can use results from the computation by Kim 
et al. [11]. They used a model 12] in which the <j> field 
was the primary order parameter undergoing phase tran- 
sition, and not a composite of the SDW order as in our 
theory above. In such a primary <fi theory, the A 2 coupling 
is strongly relevant at the decoupled fixed point [12[ (in 
contrast to the nearly marginal A2 in our theory), and a 
RG analysis is necessary to find a fixed point at non-zero 
A2. No such fixed point was found in the 3 — d expansion 
of Ref . [12I for any number of fermion flavors, while Kim et 
al. [Uj noted that a second order transition was present 
d = 2 in the limit of an infinite number of fermion fla- 
vors. Regardless of the status of this possible fixed point 
for the physical case, for our purposes here, we note that 
Kim et aVs computation of the fermion spectrum dif- 
fered from Eq. (|12p only in their assumed form for G^, 
with Eq. (dJ) replaced by G', ~ [{uj 2 + v 2 F p 2 x + v 2 A p 2 y ) 1 / 2 + 
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(uj 2 + Vpp 2 + u^p^) 1 / 2 ] -1 , the one-loop contribution of 
fermion fluctuations to the (j) propagator. Such contri- 
butions are also present in our theory, as an order A 2 , 
contribution to G^ by Gl 1 — > Gl 1 + A|G^ _1 ; given that 
A2 is nearly marginal at the decoupled fixed point, it is 
reasonable that this correction, which is formally higher 
order in A2, should be included in Eq. (fl~2|) . 

With G<p renormalized with fermion loop contributions 
as above, our results for the fermion spectral function 
have the same qualitative form as those of Kim et al. 
11]. However, it is important to note that the underly- 
ing quantum critical point, and the resulting justification 
of the computation, are very different from theirs. In par- 
ticular, an important feature of our critical point is that 
there are no restrictions on the renormalized values of 
the fermion velocities, vf and i>a, and there is no diffi- 
culty in them acquiring ratios as large as Vp/v/s, ~ 20 ob- 



served in some cuprates. As shown by Kim et al. [ll|, for 
large vp /ua > the fermion excitations are strongly damped 
by the 4> fluctuations, except in narrow arc- like regions 
around the nodal points; such "Fermi arcs" have been 
noted in a variety of experiments ljj, 2^, 21 1 . 

We can also use our results here for the scaling dimen- 
sions of the quadratic perturbations at the 0(4)® 0(3) 
fixed point to analyse the influence of quenched disorder. 
Unless magnetic impurities are explicitly introduced, the 
quenched disorder in the cuprates is spin rotation invari- 
ant, and we assume so in our discussion here. Such disor- 
der can only couple to operators, Od, which transforms 
as (D4, 13). At quadratic order in tp a i, the only allowed 
values are D = 1,9, and the scaling dimension of the 
associated couplings were quoted in Eqs. (|8l9p . After an 
average over disorder in a replica analysis, the perturba- 
tion to £J can be written as 



7d 



(13) 



where £, m are replica indices. A computation of scal- 
ing dimensions now shows that the maximum value of 
dim[7D] is for D = 9, and so the most relevant disorder- 
induced coupling is 79 with 



dim[ 79 ] =2dim[A 2 ] + 2, 



(14) 



which is clearly strongly relevant. One of the compo- 
nents of Og is (j>, the Ising nematic order, while the re- 
maining 8 components correspond to the complex charge 
density wave order parameters Q 2 , $2, $1 • $2, $1 ■ $2 
at the wavevectors 2Ki, 2K2, Ki + K2, Ki — K2 re- 
spectively. Quenched disorder will therefore induce these 
orders in glassy configurations, as was studied numeri- 
cally in Refe.EaHi. 

This paper has shown that the well-established SDW 
order of the cuprate superconductor series [l|, 0, H, @, 0] 
has precursor quantum fluctuations which can naturally 



explain key features 18, ljj, 2^, 21 1 of the spectrum of 
single electron excitations. We showed that an Ising ne- 
matic order, associated with fluctuations which reduce 
the square lattice symmetry to a rectangular symmetry 
[lot , was uniquely selected as a composite of the SDW 
order which couples most efficiently to the low energy 
single-electron excitations. 
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